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Thermal regime of the baryons behind shock waves arising in the process of virialization of dark matter halos is governed 
at cetrain conditions by radiation of HD lines. A small fraction of the shocked gas can cool down to the temperature of the 
cosmic microwave background (CMB). We estimate an upper limit for this fraction: at z = 10 it increases sharply from 
about (?y ~ 10"^ for dark halos of M = 5 x IO'^Mq to ~ 0.1 for halos with M = 10* Mq. Further increase of the halo 
mass does not lead however to a significant growth of qr - the asymptotic value for M S> 10* M© is of 0.3: We estimate 
star formation rate associated with such shock waves, and show that they can provide a small but not negligible fraction 
of the star formation. We argue that extremely metal-poor low-mass stars in the Milky Way may have been formed from 
primordial gas behind such shocks. 
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Formation of the first stars is determined by energy losses 
in ro-vibrational lines of molecular hydrogen H2 and its 
deuterated analogue HD (Lepp & Shull 1983; Shchekinov 
1986; Puy et al. 1993; Palla, Galh & Silk 1995; Galli & 
Palla 1998, 2002; Stancil, Lepp & Dalgarno 1998; Tegmark 
et al. 1997; Puy & Signore 1997, 1998; Uehara & Inutsuka 
2000; Flower 2002; Nakamura & Umemura 2002; Machida 
et al. 2005, Nagakura & Omukai 2005). In turn, the amount 
of H2 and HD and their ability to cool gas greatly depend 
on dynamical and thermal regime of the gas. In particu- 
lar, shock waves strongly enhance the rate of conversion of 
atomic hydrogen to its molecular form (Shchekinov & En- 
tel 1983; Suchkov, Shchekinov & Edelman 1983; Shapiro 
& Kang 1987; Kang & Shapiro 1992; Shchekinov 1991; 
Ferrara 1998; Yamada & Nishi 1998; Uehara & Inutsuka 
2000; Cen 2005; Machida et al. 2005; Shchekinov & Va- 
siliev 2005; Johnson & Bromm 2005). On the other hand, 
the first stars have formed from the gas processed by the 
shock waves unavoidable in the process of viriahzation of 
dark matter halos (Shapiro 1993; Haiman, Thoul & Loeb 
1996; Tegmark et al. 1997; Abel, Bryan & Norman 2000, 
2002). Therefore possible enhancement of H2 and HD in 
these conditions can have important consequences for char- 
acteristics of the first stars (Oh & Haiman 2002; Johnson & 
Bromm 2005; Shchekinov & Vasiliev 2005). 

When dark matter halos merge, shock waves form and 
compress the baryons (see discussion in (Barkana & Loeb 
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2001; Ciardi & Ferrara 2004)). At sufficiently large veloc- 
ities of colliding flows (v > 8 km s~^) fractional ioniza- 
tion in shocked gas increases above the frozen cosmologi- 
cal value, and speeds up chemical kinetics of H2 molecules. 
Collisions with velocities above v > 8.6[(1 + z)/20]~-'^/^ 
km s^^ may lead to a rapid formation of HD and an efficient 
cooling down to the minimum temperature T — Tcmb = 
2.7(1+2) (Vasiliev & Shchekinov 2005; Johnson & Bromm 
2005; Shchekinov & Vasiliev 2005). It is worth noting that 
HD can also play major role in central regions of collapsing 
dark halos, where a high-density environment favours tran- 
sition of D into molecular form, so that even in low-mass 
halos HD cooling can dominate (Ripamonti 2007). 

Explicit answer of how large is the fraction of cold gas 
in merging halos can be found only in hydrodynamic simu- 
lations. Very recently (Maio et al. 2007) performed 3D SPH 
simulation of structure formation with accounting of HD 
chemistry. They found that incorporation of effects from 
HD cooling increases dumpiness of gas, i.e. clouds are den- 
ser and more compact with respect to the case when only 
H2 cooling is considered. However, full 3D simulations are 
time consuming, and always are made within a particular 
realization of a random hydrodynamic field. They therefore 
represent only very limited regions in the space of possible 
random hydrodynamic fields of a given spectrum. This cir- 
cumstance restricts the final thermodynamic state of bary- 
ons within the domain corresponding to the chosen intere- 
lations between the amplitudes of different wave modes. As 
a result, estimates of the fraction of the cold baryons able 
to form stars are biased by such limitations. At such cir- 
cumstances ID hydrodynamical simulations of chemistry 
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and thermodynamics of shocked primordial gas can play an 
auxiliary role for understanding of what in principle can be 
expected in the conditions preceeding formation of the first 
stars. In this paper we use ID planar computations to model 
chemical and thermal regime of gas behind the shock waves 
after a head-on collision of two clouds of equal sizes. This 
approach is justified by the fact that in supersonic cloud 
collisions the rarefaction motion transverse to the symme- 
try axis takes longer than the collision time t^ = 3R/2vc 
(Gilden 1984). On the other hand, the estimates of the frac- 
tion of baryons cooled down by HD molecules we obtaine 
within this approximation are obviously upper limits. 

At the stages when the halos are close to merge their 
velocities are of the order of the virial value coiTesponding 
to the total mass of the halos. We therefore connect the col- 
Usional velocities of the clouds to the mass of the halo M 
formed in this collision 



.;e = (3Va„Pc)^/^Gi/2MV3(i + ;,)i/2 



(1) 



ilm is the matter closure parameter, pc ~ ^H'q/SttG, the 
critical density. 

In Section 2 we describe dominant thermal and chemi- 
cal processes in shocked gas; in Section 3 we discuss con- 
sequences of HD cooling for formation of the first stars; 
summary is given in Section 4. 

Throughout the paper we assume a ACDM cosmology 
with the parameters (JIq, 51a, f^m, f^h, /i) = (1-0, 0.7, 0.3, 
0.045, 0.7) and deuterium abundance 2.6 x 10^^, consistent 
with the most recent measurements (Spergel et al. 2007). 

2 Chemistry and thermal regime behind the 
shock 

In the center of mass of the colliding baryon components a 
discontinuity forms at the symmetry plane, and two shock 
waves begin to move outwards. We assume that during the 
merging collisionless dark matter components occupy con- 
siderably bigger volume and neglect gravitational forces on 
baryons. Therefore we describe propagation of the shock by 
single-fluid hydrodynamic equations. Energy equation in- 
cludes radiative losses typical in primordial plasma: Comp- 
ton cooling, recombination and bremsstrahlung radiation, 
colUsional excitation of HI (Cen 1992) H2 (GalU & Palla 
1998) and HD (Flower 2000; Lipovka, Niinez-Lopez & Avi- 
la-Reese 2005). More accurately The H2 cooling function 
have been recalculated by Le Bourlot, Pineau des Forets & 
Flower (1999), which which differs considerably (by fac- 
tor of 2) from the Le Bourlot, Pineau des Forets & Flower 
function only at T > 4000 K (Glover 2005). 

Chemical and ionization composition include a standard 
set of species: H, H+, e, H", He, He+, He++, H2, H+, D, 
D"'", D~, HD, HD+; The coiTesponding rates are taken from 
(Gain & Palla 1998;Stancilet al. 1998); the shock wave was 
computed on one collision time tc- H2 dissociation rate by 
atomic hydrogen is taken from Mac Low & Shull (1986) - 



although this reaction is important only at T > 6000 K, it 
restricts the production of H2 molecules by factor of 1 .5. 

For our simulations we use a ID planar Lagragian sche- 
me similar to that described by Thoul & Weinberg (1995). 
As a standard resolution we have used 700 zones over the 
computational region, and found a resonable convergence: 
several test computations with 1500 zones showed signif- 
icant deviations only in the very central zones of 1% in 
mass, while in the rest the deviation is normally less than 
10 %. We assume a "top-hat" initial baryonic distribution 
in colliding halos with the density equal to the virialized 
value 187r^r2bpo(l + z)'^, while temperature is taken close 
to the cosmic microwave background (CMB) temperature 
Tf, = l.lTcMB- This corresponds to a simplified picture 
when merging halos are already compressed to their virial 
radii, but not yet virialized, i.e. systematic large scale mo- 
tions did not converged into thermal energy. This assump- 
tion allows to better understand the role which shock com- 
pression of the baryons plays in their ionization and chem- 
ical state in the course of virialization. The fractional ion- 
ization X, and the abundances of H2 and HD molecules be- 
fore the shock are taken equal to their background values 
X = 10-^ /(H2) = 10-6 and /(HD) = IQ-^". 

Fig. 1 shows typical distributions of the temperature, 
the abundances of H2 and HD and their relative contribu- 
tions to the total cooling behind the shock front at t = 
(0.2, 0.6, l)tc, for the halos merged at z = 20 with the 
velocity Vc = 20 km s^^, corresponding to the total mass 
M = 1.4 X IO^Mq. Three qualitatively different cooling 
regimes can be distinguished in the temperature profiles: 
in the high temperature range (1500 < T < 7000 K) ex- 
citation of ro-vibrational levels of H2 dominates, while in 
the intermediate range (200 < T < 1500 K) only rota- 
tional Unes contribute to the cooling; in the lowest temper- 
ature range (T < 200 K) cooling from H2 molecules ex- 
hausts and only HD rotations support cooling - it is seen in 
the lower panel from comparison of the relative contribu- 
tions of H2 and HD cooling. At later times, t ^ tc a non- 
monotonous behaviour of temperature is due to the fact that 
at these stages gas enters the shock front with lower veloci- 
ties, so that the post-shock temperature becomes below the 
limit when collisional dissociation of H2 molecules is im- 
portant. As a result, all molecular hydrogen formed behind 
the front survives and stimulates rapid cooling. 

In this particular case a small fraction (q ~ 0. 1 by mass) 
of the shocked baryons close to the symmetry plane cools to 
the minimum possible value T ~ Tcmb ~ 2.7(1 + z) due 
to cooling in rotational lines of HD. In general, the frac- 
tion of compressed baryons cooled down to a given tem- 
perature level depends on the relative velocities of the col- 
liding clouds Vc'- the larger the collisional velocity Vc, the 
stronger the gas compression after cooling, and the higher 
the contribution from HD cooling. Fig. 2 shows the fraction 
of baryons qt {M) contained in several temperature ranges 
versus the halo mass: T < 200 K, T < 150 K, T < 100 
K, and at T ~ Tcmb after one crossing time t = tc- In 
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Fig. 1 Upper panel: profiles of temperature (solid) and den- 
sity (dash); middle panel: relative concentration of H2 and HD 
molecules; lower panel: their relative contribution to the total cool- 
ing (H2 - solid and HD - dashed), for baryons in two colliding 
halos with the total mass M = 1.4 x lO'^A/0 at 0.2tc, OMc, tc\ 
halos merged at z = 20. 



the temperature range T < 200 K where cooling in H2 
lines is dominant, qT{M) increases from 0.1 for the halo 
mass l{fMQ to the asymptotical (at 7\/ > 2 x l(f Mq) 
value 0.6. Compressed baryons can have temperature be- 
low 150 K only due to contribution from HD cooling. In the 
lower temperature ranges (T < 100 K, and at T = Tcmb) 
qriM) is a very sharp function of the halo mass: for in- 
stance, at redshift z = 20 (Fig. 2a) a two-fold increase of 
the mass from lO^M© to 2 x l{f Mq results in a two-order 
of magnitude increase of Qtcmb from 3 x 10"'^ to 0.3. At 
higher masses the dependence flattens and asymptotically in 
the limit 7\/ » 2 x IO^Mq approaches 0.4. At lower red- 
shifts gas density decreases and the radiation cooling time 
becomes longer As a result, qT{M) shifts towards larger 
halo masses, approximately as (1 + z)^^/^: Fig. 2b shows 
qriM) for collisions occured at z = 10. 
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Fig. 2 Fraction of baryons cooled below temperature T < 200 
K, r < 150 K, T < 100 K, and T = Tcmb from top to bottom 
at 2 = 20 (left), and z = 10 (right). 



3 Star formation 

Baryons cooled below T < 200 K are thought to be able 
to fragment and initiate star formation. In gas layers com- 
pressed by shock waves gravitational instability of the cold 
gas occurs naturally when the thickness of the layer equals 
the Jeans length. In order to estimate the range of masses 
expected to form through the gravitational instability we ap- 
ply (Gilden 1984) criterion for the shock-compressed gas, 
which implies that: /) the characteristic growth time is short- 
er than the collision time, and //) the critical wavelength is 
shorter than the initial size of the clouds. The correspond- 
ing critical mass A/cr depends on the average temperature 
and density in the layer: when the halos merge with low rel- 
ative velocities (corresponding to lower halo masses), only 
a small fraction of the compressed baryon mass can cool 
down to sufficiently low temperatures to form an unstable 
layer Mergings with higher relative velocities increase the 
fraction of gravitationally unstable baryons. Fig. 3 shows the 
redshift dependence of the halo masses with a given fraction 
of the compressed baryons unstable in Gilden sense. Each 
line is marked with symbols corresponding to a given frac- 
tion of the baryon mass qj unstable against fragmentation: 
for instance, when halos with masses corresponding to the 
upper line M = 2.5 x 10'^[(1 + z)/20]^'^-'^Mq merge, half 
of their mass becomes compressed in a cold layer with tem- 
perature r < 100 — 150 K unstable in Gilden sense. Co- 
incidently, the thermal evolution is governed by HD cool- 
ing, when the temperature falls below 150 K. At the lat- 
est stages the unstable layers are dominated by HD cool- 
ing, so that fragments formed in these conditions can reach 
the minimum possible temperature c± Tcmb- The corre- 
sponding Jeans mass in the unstable layer is Mj < 2.3 x 
lO^M0w[o^[(l + z)/20]i/2, which is considerably smaller 
than the baryonic mass of this layer; here wio = t'c/10 
km s~^. With accounting (1) the Jeans mass in the unsta- 
ble layer is Mj < 1.3 x IO-^A/^^/^Mq. This means that 
when halos with masses A/ > lO'*A'/0 merge more than 
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MF{M) dM, 



(2) 



Fig. 3 Lines with simbols depict halo masses vs redshift, where 
the fraction of baryon mass unstable in Gilden sense is at given 
level: 1%, 7%, 25%, 50% from bottom to top; straight thin solid 
line corresponds to a 3a peak mass, dashed line shows the min- 
imum mass obtained by (Tegmark et al. 1997), thick solid line 
shows the minimum mass from (Shchekinov & Vasiliev 2006). 



one cold and dense clouds in the unstable layer can form 
and give rise to formation of stars. 

Whether such fragments are the protostellar condensa- 
tions evolving further in a single massive star, or they are 
the nestles, where in the process of sequential (hierarchi- 
cal, Hoyle 1953) or of a one-step (Nakamura & Umemura 
1999) fragmentation a cluster of less massive stars is form- 
ed, depends sensitively on details of gravitational contrac- 
tion and thermal evolution of the fragments (see discussion 
in (Coppi, Bromm & Larson 2001; Glover 2004). However, 
both the mass of a central protostellar core in the former 
case, and the minimum mass of a protostellar condensa- 
tion in the latter are determined by the opaqueness of the 
contracting gas. Vasiliev & Shchekinov (2005) estimate this 
mass as Mj ~ 10~'^(1 + z)^^^Mq, what gives for z — 
10 - 20 a relatively low mass limit: Mj ~ (0.03 - O.1)M0. 
Therefore stars formed in such cold layers with a predomi- 
nance of HD cooling on later stages, are anticipated in gen- 
eral to be less massive than those expected when thermo- 
dynamics of primordial gas is determined by H2 cooling 
(Shchekinov & Vasiliev 2005). 

The fraction of baryons in the universe able to cool be- 
low T = 150 K and to form presumably low-mass stars 
behind shock waves in mergings can be estimated as 



JV/„. 



where F{M) ~ dN /dM is the Press-Schechter mass func- 
tion, Mc is the halo mass at a given redshift presumably 
formed through merging of two M^/l halos, where the frac- 
tion of cold (T < 150 K) baryons is equal to 1%, Mmin, the 
minimum halo mass at a given redshift. We assume that the 
halos with the fraction of cold baryons less than 1 % do not 
significntly contribute to /c or the star formation rate. Here 
we introduced factor €±_ = 0.05 accounting only approx- 
imately head-on collisions. Indeed, our conclusions about 
the role of HD cooling are based on the assumption of a 
head-on collision of merging halos, and can be valid only 
in a restricted range of the impact parameter when the shear 
motion is less important than the converging flow and the 
corresponding diverging shock waves. For this condition to 
be fulfilled the characteristic time of the Kelvin-Helmholtz 
instability of the shear flow t^^ ~ R/v\\ must be longer 
than the dynamical (crossing) time td ^ 2R/v±: t-^-^ > td- 
This gives v\\/v±_ < 1/2, and as a result, only a fraction 
ex = Ail/47r ~ 0.05 of mergers where the flows are ap- 
proximately head-on. In Introduction we stressed that the 
estimate of the cold baryon fraction /c within the approxi- 
matrion of head-on collisions is an upper limit. Introducing 
the factor e^ attempts to quantify more accurately the frac- 
tion fc- However, the estimate (2) still has a meaning of an 
upper limit, particularly because /) in (2) we assumed ex- 
plicitly that the spectrum F{M) above Mc is formed thro- 
ugh merging of equal masses, and //) we do not account here 
(negative) feedback from the stars have already formed dur- 
ing the merging. 

Fig. 4 depicts fc versus redshift. For M„iin we have 
taken IQ^Mq, IQ'^Mq and 10^ A/©, however, only the last 
value seems meaningful, because for lower halo masses the- 
ir baryonic content is too small and likely can be easily re- 
moved in tidal interactions. Independent on A/min the total 
fraction of baryons in the universe able to cool below T = 
150 K and form stars approaches fc ~ O.le^ at z = 10. 
This fraction may increase if the halos are clustered as sug- 
gested by Barkana & Loeb (2004) and Cen (2005). For com- 
parison we show in Fig. 4 the fraction of baryons in the uni- 
verse cooled by H2 molecules and hydrogen atoms with the 
lower limit Mc in the integral (2) corresponding to the virial 
mass Mc = My at T^i,. = 400, lO'* K, for both cases we 
assume that in virialized halos 8% of baryons can have tem- 
perature 200 K as calculated by (Abel et al. 1998), and inte- 
grate over the halo mass spectrum from Afmin = 10'* A/0. 

Characteristic star formation time inside the formed fra- 
gments is determined by the baryon density in there igf = 
tj/e, where tj ~ ^3tt/32Gp is the baryon Jeans time in 
cold layers, e ^ 1 is the star formation efficiency. Gas den- 
sity pf in cold layers when they fragment can be found from 
the condition piVc = kpfTcMB/mH, which gives 



V 



0.02Af7/^(l 



z)^ cm ^, 



(3) 
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Fig. 4 Fraction of baryons fc in the universe cooled below 
150 K. According to Gilden criterion these baryons can give rise 
to the formation of stars. Solid curves show the fraction fc for 
Mmin = 10^ M0, 10* Mq, W^Mq (bottom to top) in eq. (2). 
Dashed curves depict two cases when the lower limit in the r.h.s. of 
(2) is equal to the virial mass Mc = M„ with T^ir = 400, lO" K 
(top to bottom); the lower limit Mmin in the denominator of the 
r.h.s. of (2) is taken for these two curves to be W^Mq. 



and 



ij = 5 X 10* A/7 ^^^(1 



z)-3/Vr, 



(4) 



M,^--Lee^ / dMiF(Afi) 



For star formation efficiency e > 0.03 the charactersitsic 
time isf remains shorter that the Hubble time for the halo 
masses Mh > 10^ Mq. This means that star formation rate 
in merging halos is determined by the longest time - the 
charactersitic time between subsequent mergings. In these 
conditions star formation rate is proportional to the merger 
rate of the halos with the total mass above the critical value 
Mc{z) (Barkana & Loeb 2000; Santos, Bromm & Kamion- 
kowski 2002) 

-A/i 

dM.,,(M..)A.f.^, (5) 

where P = P{Mi,M2,z) is the probability that a halo 
with mass A/i merges to a halo of mass M2 > Mi at 
redshift z (Lacey & Cole 1993); we explicitly assume here 
that only fraction of baryons qf{M, z) cooled after merging 
below 200 K is able to form stars. Therefore, if we sub- 
stitute here the fraction qT<i5o{M) of baryons cooled to 
T < 150 K, equation (5) will describe the contribution to 
the total star formation rate from the merging halos where 
thermodynamics is governed by HD cooling, and where low 



Ma 



mass stars may form. It is obvious that mergings of halos of 
masses A/min ^ A/c with the halo of critical (or overcriti- 
cal) mass, involve too small baryon mass fraction into suf- 
ficiently strong compression where HD molecules can cool 
gas. Moreover, in this limit the compressed region deviates 
significantly from planar geometry, and Gilden criterion is 
not applicable. Therefore, in our estimates of star formation 
rate we assume for Af,ui,i two values: M^in = 0.5 A/^, and 
-^^min = 0.9 Mc with Ale defined as above as the mass 
where the fraction of cold baryons {T < 150 K) after merg- 
ing is equal to 1%; at equal densities of merging halos their 
sizes differ in the first case by ~ 20%, while in the second 
only by 4%. Therefore, in the estimates of the star formation 
rate with a predominance of HD in thermodynamics, we 
neglect contribution from mergers of masses significantly 
smaller than the critical mass. On the other hand, it is seen 
from Fig. 3 that halos with ^ 2A/c have already ~ 50% of 
their mass colder than T < 150 K and unstable in Gilden 
sense. This means that the upper limit 2Afc in the second 
integral of (5) counts practically all halos whose star forma- 
tion is regulated by HD cooling. With this proviso, the two 
cases; A/,nin = 0.5 Mc and A/,nin = 0.9 A/c, are shown in 
Fig. 5 by two solid lines - the region between them can be 
a reasonable estimate of the star formation rate governed by 
HD cooling. 

For comparison in Fig. 5 we add two lines for which 
A/c in (5) is replaced by the virial mass My with Tyir — 
400, 10^ K (see, Barkana & Loeb 2000). It is obvious, that 
the number of mergers where HD cooling dominates, is only 
a small (although not negligible) fraction of all mergers; 
at z between 10 and 16 star formation governed by HD 
varies from 15 to 30 % of the one connected with 10^ K 
halo mergers (the lower dashed line), and 0.5 to 10 % of 
star formation in 400 K mergers (the upper dashed line); 
in Fig. 5 for all low-mass mergers dominated by H2 cool- 
ing 9/ (A/, z) = (7/.H2 — 0.08 is assumed following (Abel 
et al. 1998). At earlier stages, z = 18 — 20, mergers with 
a predominance of HD cooling contribute less than 0.5% 
compared to the 400 K mergers. From this point of view 
one can think that in numerical simulations the regions with 
HD cooling are apparently missed. 

It is therefore seen that a small fraction of baryons in 
mergers cools down to the lowest possible temperature T ~ 
TcMB, and can give rise to formation of the first genera- 
tion stars of low masses - lower than the masses formed 
under the conditions when H2 controls thermal evolution 
of baryons. This fraction increases with the total mass of 
merging halos, and thus in massive galaxies the popula- 
tion of low-mass stars of the first generation can be con- 
siderable. One should stress though that massive galaxies, 
formed in the hierarchical scenario through mergers, are 
quite expected to have been already experienced star forma- 
tion episodes, with the insterstellar gas polluted by metals. 
However, it remains unclear whether the metals can become 
well mixed in a galaxy before it absorbs a new halo in next 
merger event. Simple arguments suggest the opposite. In- 
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Fig. 5 Cosmic star formation rate: dasiied lines correspond to 
the halos with T^ir = 400, 10* K (top to bottom); solid - the 
halos with mass Mcr ~ Mc, Mmin = 0.5Mc (upper curve) and 
Afmin = 0.9Mc (lower curve); in all cases e = 0.1. 



deed, the characteristic mixing time for the whole galaxy 
can be estimated as tmx ~ {Sej)R/cs, where (S^j) > 1 is 
the mean density contrast between SNe ejecta and diffuse 
interstellar gas, R is the galaxy radius, Cg, the sound speed 
in diffuse gas; note that R/cs ^ tc- The characteristic time 
scale for the halo mass growth tm ~ [M2cPP{Mi, A/2, z)/ 
dAhdt]-^. For the halos M2 ~ 10'^ Mq and Mi - O.9M2 
at z = 20, t„i is about ^ 8 x 10^^ s, which is comparable 
to the collision time tc = 3R/2vc, and is therefore < tmx- 
Although some of the absorbed low mass halos can have 
been also experienced star formation episodes before being 
merged, and thus can be already metal enriched (Scanna- 
pieco, Schneider & Ferrara 2003), however a non-negligible 
fraction of them may have pristine composition due to a 
slow mixing. 

The existence of low mass Pop III stars is suspected 
from the observational point of view: recently discovered 
extremely metal-poor and low-mass stars, as for instance a 
0.8 Mq star HE 0107-5240 (Christlieb et al. 2002, 2004) 
with [Fe/H]= —5.3, may be the first-generation stars. The 
question of whether stars with [Fe/H]< —5 are indeed Pop- 
ulation III stars is still under discussion, partly because of 
overabundant carbon and nitrogen: [C/Fe]= 4, [N/Fe]~ 2.3 
in HE 0107-5240 (Bessel, Christlieb & Gustafsson 2004). 
From this point of view HE 0107-5240 can be a Population 
II star formed in an already enriched interstellar gas (Umeda 
& Nomoto 2003). However, the possibility that this star has 
been formed of pristine gas cannot be excluded: Shigeyama 
Tsujimoto & Yoshii (2003) conclude that HE 0107-5240 is 
a Pop III star with the surface polluted by accreted interstel- 
lar gas already enriched with metals. In this scenario over- 
abundant C and N in the envelope can be produced during 



the core helium flash as suggested by (Weiss et al. 2000; 
Schlattl et al. 2002). The estimates presented in this paper 
show that a small fraction of baryons processed by shocks in 
merger events may indeed become extremely cold and form 
low mass stars of the first generation. 

4 Conclusion 

In this paper we have shown that 

- A small fraction of baryons in merging halos can cool 
down to very low temperatures close to the temperature 
of the cosmic microwave background. This fraction in- 
creases with the halo mass, and can reach ~ 0.2 for 
masses Af > 2 x 10^ Mq at z = 20 and for M > 
10^ A/0 at 2 = 10; averaged over the halo mass spec- 
trum this fraction varies from 10^"^ at z = 20 to 0.1 at 
z = 10. 

- Such cold gas is unstable against gravitational fragmen- 
tation, with the mass of the primary fragments decreas- 
ing for mergers of higher masses M: Mj < 1.3 x 10^ 
AI~^^^Mq. Masses of protostars formed in gas cooled 
by HD molecules are in general lower than those formed 
in conditions when H2 cooling dominates. 

- The contribution to the cosmic star formation rate of the 
mergers with a predominance of HD cooling, and there- 
fore with presumably low mass first stellar objects, in- 
creases from less than 0.5% at redshift z = 18 — 20 to 
10-30 % at z = 10. Extremally metal-poor low mass 
stars in the Milky Way may have been formed in merg- 
ers dominated by HD cooling. 
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